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Abstract. Let M° be a complete noncompact manifold and g an asymptot- 
ically conic manifold on M° , in the sense that M° compactifies to a manifold 
with boundary M in such a way that g becomes a scattering metric on M. A 
special case that we focus on is that of asymptotically Euclidean manifolds, 
where the induced metric at infinity is equal to the standard metric on S n_1 ; 
such manifolds have an end that can be identified with R n \ B(R, 0) in such 
a way that the metric is asymptotic in a precise sense to the fiat Euclidean 
metric. We analyze the asymptotics of the resolvent kernel (P + k 2 )~ 1 where 
P = A 9 + V is the sum of the positive Laplacian associated to g and a real 
potential function V that is smooth on M and vanishes to third order at the 
boundary (i.e. decays to third order at infinity on M°). We show that on a 
blown up version of M 2 X [0, fco] the resolvent kernel is conormal to the lifted 
diagonal and polyhomogeneous at the boundary, and we are able to identify 
explicitly the leading behaviour of the kernel at each boundary hypersurface. 
Using this we show that the Riesz transform of P is bounded on L P (M°) for 
1 < p < n, and that this range is optimal if V ^ or if M has more than one 
end. The result with V ^ is new even when M° = R n , g is the Euclidean 
metric and V is compactly supported. When V = with one end, the range 
of p becomes 1 < p < p max where p max > n depends explicitly on the first 
non-zero eigenvalue of the Laplacian on the boundary dM. 

Our results hold for all dimensions > 3 under the assumption that P has 
neither zero modes nor a zero-resonance. In a follow-up paper we shall analyze 
the same situation in the presence of zero modes and zero-resonances. 



We study the resolvent (P+fc 2 ) -1 as k J, for Schrodinger operators P = A g + V 
on asymptotically conic manifolds. Let M° be a smooth manifold with one or more 
ends, each diffeomorphic to Si x (0, oo) where Si is a smooth compact manifold 
of dimension n — 1. Write S for the disjoint union of the Si. We suppose M° 
admits a compactification M which is smooth with boundary dM — S and let g 
be an asymptotically conic metric (also called scattering metric) on M, that is a 
complete metric on M such that in a collar neighbourhood [0, e) x x dM near dM, 



where a: is a smooth function that defines the boundary dM (i.e. dM = {x = 0} 
and dx does not vanish on dM) and h(x) is a smooth family of metrics on 5". We 
let V be a real potential such that 



A particular case which will be of special interest in this work is the case where 
S = 5 f,i_1 is the sphere and ho '■= h(0) = d8 2 is its canonical metric. Then we say 
that g is asymptotically Euclidean since the Euclidean metric in polar coordinates 
reads dr 2 + r 2 d9 2 on R+ x S 11 ^ 1 and can be put under the form dx 2 /x 4 + d9 2 /x 2 
by setting x := 1/r. For a general scattering metric taking the form we can 



1. Introduction 





(1.2) 



VeC°°(M), V = 0(x 3 ) as x -> 0. 
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view r — 1/x as a, generalized 'radial coordinate', as the distance to any fixed point 
of M is given by r + 0(1) as r — » oo. 

Let A g be the positive Laplacian on M with respect to g. Then the Schrodinger 
operator 

P = A g + V 

is self-adjoint on L 2 (M,dg) and its spectrum is cr(P) = [0,oo) U a pp (P) where 
a pp(P) = ^ ' ' ' ^ — ^tv! i s a finite set of negative eigenvalues (by convention 

ki > 0). Note that can be an L 2 -eigenvalue but the half- line (0, oo) is only con- 
tinuous spectrum. 

The resolvent R(k) = (P + k 2 )^ 1 is well defined as a bounded operator on 
L 2 (M, dg) for k 6 (0, fci) but fails to be bounded or defined at k = 0. We describe 
in this work the uniform behaviour near k = in terms of distributional kernel of 
R(k; m, m') on the manifold with corners X = [0, fco]/c x M x M for small ko < k±. 
More precisely, we follow Melrose's program [15] . pass to a desingularized version 
of X, that we call M 2 sc (see figure E2T]), obtained by blowing up a certain number 
of corners of X, and lift the kernel of R(k) to M? sc . Our result is that R(k) is 
conormal to the lifted diagonal of M 2 sc and polyhomogeneous conormal at the 
boundary. Thus, when lifted to M 2 sc the kernel has asymptotic expansions at 
all orders at each face. In particular, the expansion at the face corresponding to 
M 2 x {k — 0} is simply the Laurent expansion of R(k) at k — 0. 

This Laurent expansion is quite important if one wants to understand the large 
time asymptotics of Schrodinger propagators, heat operators or wave operators. 
It has been studied by Jensen-Kato [5] for (M°,g) = (R 3 ,g euc i) and extended by 
Murata [H] for M° = M" and P = p(D) + V a real constant coefficient elliptic 
differential operator with V a compact operator in weighted Sobolev spaces. Then 
more recently X.P. Wang [22] studied the case of the Laplacian plus a potential 
which is homogeneous of degree —2 on an exact conic manifold, perturbed by a 
shorter range potential. He showed there is an asymptotic expansion of R(k) at 
k = with a number of terms that depend on the decay rate of the perturbation. In 
all these approaches, an exact formula for the free model operator is used to derive 
the Laurent expansion and compute its leading coefficients, which are expressed in 
terms of the model resolvent. It is not clear how to generalize this approach to 
cases with no free model. 

Our approach is to constuct the resolvent kernel through a parametrix, guessing 
what must be the Taylor coefficients at each face of M 2 sc and in particular at k = 0. 
This, in turn, relates to the analysis of a Laplacian type operator on a manifold 
(M,gb) conformal to (M,g), 

dx 2 

g b := x 2 g = — + h(x), 
x z 

that is, a manifold with asymptotically cylindrical ends. We use the b-calculus of 
[TT] to solve the model problems at k = and, as a corollary, we describe the first 
coefficients of the Laurent expansion of R(k) at k = 0. Our results are phrased in 
terms of a calculus of pseudodifferential operators \j/™^ abf o ,azf,aBC -' !£ ^/; f2^ 2 ) act- 
ing on half-densities on M, where m is the pseudodifferential order, (<ibf , a z f, a sc ) 
determines the growth of the Schwartz kernel of the operator at the three boundary 
hypersurfaces meeting the diagonal (see Figure 12". 2. lj) and £ is an 'index family' list- 
ing the allowable terms in the asymptotic expansion of the kernel at each boundary 
hypersurface of M 2 SC . 
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Theorem 1.1. Suppose that n > 3 and that P has neither zero-modes nor a zero- 
resonanc^. Then 

(1.3) R(k) G $-2,(-2,0,0),S (M! fjl/2 ) 

where 3? is an index set. This index set is given explicitly in (|3.10[) w/ien M is 
asymptotically Euclidean and lower bounds are given in (|4.3[) for the general case. 
More informally, R(k) is conormal to the diagonal of M| sc and polyhomogeneous 
at the boundary, uniformly down to k = 0. 

We use this result to analyze the boundedness of the so-called Riesz transform 
as an operator from L p (M,dg) to L p (M,T*M,dg). If n> denotes the spectral 
projection for the operator P onto the positive spectrum (0, oo), and P> = Pon> 
then we define the Riesz transform T by 

T :=do(P > )-2. 

Here d is exterior differentiation from functions to 1-forms. 

The Riesz transform for the Laplacian on a complete non-compact manifold is 
clearly bounded on L 2 but its L p boundedness is in general quite involved and can 
be interpreted as a way of comparing two definitions of Sobolev spaces. We refer the 
reader to the paper [3] for general results and its section 1.3 for the state of the art 
about Riesz transform LP boundedness. Here we mention that for the Laplacian, 
T G £j(L p ) if p £ (1,2] in a quite general setting (see Coulhon-Duong 0); in other 
words low ranges of p seem less sensitive to the geometry. For p > 2, the question 
is more delicate and the first cases to look at are manifolds which have a Euclidean 
structure near infinity. Using the formula 

p> 2 = - / i?(/c)n> dk 

n Jo 

Carron-Coulhon-Hassell [3] analyzed the LP boundedness properties of T in this 
setting (again V = 0) through a careful analysis of the resolvent. They showed 
that T 6 L(L P ) if and only if p g (1, oo) when the manifold has one end, whereas 
p G (1, ri) when there is more than one end. The unboundedness for p > n had 
been shown earlier by Coulhon-Duong [6]. 

As for the case V ^ 0, some work has been done on LP estimates, especially 
for Schrodinger operators in M. n , but in most cases, it seems that the potential is 
supposed non-negative so heat kernel techniques are effective. We note for instance 
the boundedness of Riesz transform on LP for p G (1,2] if V > and V G L 1 1 oc (R") 
(see E.M. Ouhabaz [19]), which means again a quite strong stability for low range 
of p. A similar result is obtained by Coulhon-Duong [6 for the Laplacian on a class 
of manifolds. P. Auscher and A. Ben Ali [2] improved a result of Z. Shen [SO] on 
higher ranges of p, this is when the potential is non- negative and satisfy the q reverse 
Holder inequality for some q G (1, +oo], indeed they prove in this setting that there 
is e > such that T G L(L P ) if p G (1, q+e) (Z. Shen's restriction was n/2 < q < n). 
We remark that powers \z\~ a satisfy the reverse Holder inequality if a G (— oo, n/q). 
It is interesting that their hypothesis rules out potentials V decaying too fast near 
infinity, while our hypothesis by contrast requires sufficiently fast decay at infinity. 

Continuing the program of Carron-Coulhon-Hassell [1] (see their last section), we 
extend their results to more general manifolds and include the case with potential, 
which need not be non-negative. We do assume in this paper that P has neither 
zero modes, nor a zero-resonance. We shall see, in a sequel to this paper, that the 
presence of zero modes and zero-resonances has dramatic consequences for the Riesz 

^By a zero mode we mean an L 2 eigenfunction with eigenvalue 0, and by a zero-resonance a 
solution v to Pv = that decays to zero at infinity. Zero-resonances can only exist for n < 4. 
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transform L p boundedness. This analysis in turn should facilitate the analysis of 
the Riesz transform for differential forms on asymptotically conic manifolds, since 
then L 2 harmonic forms and resonances may appear. 

We first state our Riesz transform results in the case V = 0. 

Theorem 1.2. Let n > 3, and let P be the Laplacian on an asymptotically conic 
manifold [M, g) with one end. Let Aqm be the Laplacian on the boundary of M 
for the metric h(0) given by il.l]) , let Ai be its first non-zero eigenvalue, and let 
v\ = y/ ((n — 2)/2) 2 + Ai. If V\ < n/2, then the Riesz transform T satisfies 

T e L(L P (M),L P (M;T*M)) Kp<n(|-n) 

while if v\ > n/2, then 

T e L(L P (M), L P (M;T*M)), for all 1 < p < oo. 

Notice that n(n/2 — Vx)^ 1 , the upper threshold for p when z/i < n/2, is always 
greater than n. This result answers Problem 8.1 of [4]. It extends that of [4] for 
Euclidean ends and is closely related to (and indeed uses) that of Li [TUj, who 
proves LP boundedness in the same range for the Riesz transform on metric cones. 
One direction of these equivalences can also be obtained by the recent paper of 
Coulhon-Dungey [BJ with a different method. Next, we show 

Theorem 1.3. Let n > 3 and let P = A g + V with V satisfying il.ty) and such 
that P has no zero modes or zero-resonance. Suppose that either M has more than 
one end, or that V ^ 0. Then 

T e L(L p (M),L p (M; T*M)) <^ 1< p < n. 

Remark 1.4. This seems surprising at first, since P may be a very mild perturbation 
of the flat Laplacian on R™ (if V is, say, compactly supported and nonnegative). 
Yet it reduces the upper threshold for L p boundedness of T from oo to n. The 
lack of boundedness for the Riesz transform for p > n for the case of the connected 
sum of two copies of R" was remarked on by Coulhon and Duong [BJ, but there 
it appeared to be the global geometry and topology of M that was responsible for 
this phenomenon. However, [4] and the present paper give a different perspective; 
it is whether there is a bounded nonconstant solution to Pv = that determines 
whether the Riesz transform is unbounded for p > n. 

Remark 1.5. To obtain these Riesz transform results it is essential that we under- 
stand the behaviour of R(k) at all the boundary hypersurfaces at k = 0, and not 
just at the lift of M 2 x {0}. Neither Jensen-Kato nor Murata do this, and conse- 
quently, it does not seem possible to deduce boundedness properties of the Riesz 
transform from their results. 

Remark 1.6. It can be seen from the Hardy inequality that if M = M. n with n > 2 
and such that V^_(z) < a\z\~ 2 for a < (n/2 — l) 2 where V- is the negative part of 
V, then P = Ar« + V does not have a zero-resonance nor nontrivial L 2 kernel. In 
particular, the frequently-made assumption V- — avoids such problems. 
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2. Preliminaries 

2.1. Manifolds with cylindrical ends. As mentioned in the introduction, we 
convert the problem of inverting P + k 2 to an equivalent problem on a Riemannian 
manifold with asymptotically cylindrical ends. Let us assume that M has only one 
end E ~ M+ x dM to simplify the exposition. When the manifold is asymptotically 
Euclidean, we have coordinates z = (rt/i, . . . , ry n ) on E induced by the coordinates 
y G R" on the sphere {\y\ — 1}, and we denote by x = |z| _1 = 1/r a smooth 
function that defines the boundary of the radial compactification M. We denote by 
Vb the Lie algebra of smooth vector fields on M that are tangent to the boundary 
dM and Diffb(M) the enveloping algebra of V& over C°°(M), i.e. the set of dif- 
ferential operators in these vector fields. Such vector fields can be locally written 
as a combination ovr C°°(M) of xd x , d yi , . . . , d Vrl _ 1 where (d Vi )i form a local basis 
of the tangent space of the boundary dM. We denote by V sc the set of smooth 
vector fields which can locally be written near the boundary as a combination of 
x 2 d x , xd yi , . . . , xd Vn _ 1 and Diff sc (M) the enveloping algebra of V sc . 

We define gt := x 2 g the conformal metric to g and L 2 (M) = L 2 (M, dvol 9i) ) = 
L 2 (M). (M, <7b) is then an exact b-metric in the sense of Melrose [IT], that is 
a cylindrical metric on M . We now define the operator Pb by 

(2.1) P = x^ +1 P b x-* +1 . 

Since P is formally self-adjoint with respect to g, Pb is formally self-adjoint with 
respect to gb- A calculation gives 

/ r? — 2 \ 2 

(2.2) P b := -{xd x f + J +A 9M + W, 

with W £ xDiSb(X) a lower-order term at x = 0. 

We study P indirectly via analyzing the operator P& using the b-calculus and 
results from [llj . We now summarize the major results of the b-calculus as they 
apply to our particular operator P b . First we recall the b-Sobolev spaces Hl(M) 1 
j = 0, 1, 2, ... , consisting of those functions in L 2 (M) which are mapped into L 2 (M) 
by all order j elements of Diffb(M). The weighted b-Sobolev space x a H° b (M) is 
defined as0 the space of functions which can be written x a f with / G fl^ (M). Let 
Ao = < Ai < A2 . . . be the spectrum of the Laplacian on (dM,ho). We then 
define 

(2.3) N g := {v , I Vi = V((n - 2)/2) 2 + Ai}. 

In the case of the canonical sphere (dM — S"^ 1 , ho — dd 2 ) we have 

N d = {(n-2)/2 + l;leN Q }. 
Melrose's Relative Index Theorem (|llj. section 6.2) implies in our setting 

Theorem 2.1. The operator Pb is Fredholm as a map fromx a H 3 b (M) to x a H 3 b ' 2 (M) 
for all j > 2 and all a ivi, I = 0, 1, 2, ... . The index of Pb is equal to for 
\a\ < (n— 2)/2 and the index jumps by di, the multiplicity of the vf-eigenspace E Vl 
of Aqm + (n — 2) 2 /4, as a crosses the value ±vi (with a decreasing). 

We also need results on the regularity of solutions of the equation PbU = f, when 
/ is poly homogeneous. The following result proved in (llj is phrased in terms of 
index sets and the operation of extended union; see subsection 12.2.41 



The b-Sobolev spaces can be defined for any real j, but it is not necessary to do so for our 
purposes. 



6 



COLIN GUILLARMOU AND ANDREW HASSELL 



Theorem 2.2. Suppose that f is polyhomogeneous on M with respect to the index 
set E, that u G x a L 2 (M), and that P^u = f . Let S be the set 

S = {± Vl | 1 = 0,1,2,...} 

and for b G R let fJ-(b,j) — + k;k — 0,...j}nS. Then u is polyhomogeneous 
with respect to the index set EUF, where F is the index set 

{((±n + j),k) 1 1 e No, j e N , ±n>a, o<k< ^{±v u j) - 1}. 

When (dM,h ) = {S'^^dO 2 ), this reduces to 

{(n/2 + l,k) I G Z, n/2 + l>a, < k < N t - l}. 

where Ni is the number of elements of the form ±(n/2 — 1+j), j £ No in the interval 
(a,n/2 + l] . 

As a consequence of these two theorems, the vector space 

{/ G x^Hl(M) | P b f = 0}/{/ G x^Hl(M) | P b / - 0} 

for sufficiently small e > is finite dimensional, independent of e, and independent 
of j by elliptic regularity. By Theorem 12. 2[ elements of this space have the form 
x Vl 4>i{y) + 0(x Ul+e ) where 4>i(y) is an element of E Vl . The span of such elements 
4>i is a vector subspace G Ul c -EWji- 

Proposition 2.3. ([llj. Chapter 6) The subspaces G Vl and G- Vl of E Vl are or- 
thogonal complements with respect to the inner product on (dM,ho). 



Remark 2.4. Thus, the index jumps by di as a crosses the value v\ since the dimen- 
sion of the null space increases by dim G Vl , while the dimension of the cokernel, 
which can be regarded as the dimension of the null space of Pb on x~ l H^(M), 
increases by dim G- Vl . 

2.2. The blown-up space M| sc . 

2.2.1. Definition of M% . We want to analyze the distributional Schwartz kernel 
of 

R{k) := (P + k 2 )- 1 
as k — > and k > 0. The model case with X = M 3 , g = g cuc \ and V = gives 

1 e -k\z-z'\ 



R(k;z,z') = 



47r \z-z'\ 

this shows that the different asymptotic behaviours of this kernel on [0, l]fc x M x M 
leads to consider it on a blow-up space M| sc of [0, 1] x M x M as described in a 
note of Melrose-Sa Barreto [T5]. For instance R(k; z, z') goes from decreasing at all 
orders (i.e. 0(2;°°)) as x = — > when k > 0, z ^ z' to decreasing at order 1 
(i.e. 0(x)) when k = 0. On the other hand, in the asymptotic regime (x, k) — * 0, 
z' fixed in M, the kernel is a smooth function of the variables z', z/\z\, x and k/x. 
This suggests that 'polar coordinates' in (x, k) are more adapted to the analysis 
of R(k,z,z'). We realize these and other singular coordinates geometrically by 
working on the Melrose-Sa Barreto space in which they appear as smooth functions. 

We use the unpublished note [H] and refer to [T7], [TB] and [TT] for information 
and notation about manifold with corners, conormal distributions, blow-ups and 
so on. Here we recall that a submanifold Y of a manifold with corners X is said 
to be a product-type submanifold or p-submanifold if for every y G Y there exist 
local coordinates X\, . . . , x r , y' , y" with Xi > local boundary defining functions 
and y' € M s ', y" G R s " , r + s' + s" = dimX so that Y is locally defined by 
{xi = ■ ■ ■ = x q = 0, y' = 0} (q < r). We denote by [X; Y\, . . . , Y/v] the iterated real 
blow-up of X around N submanifolds Yi if Y\ is a p-submanifold, the lift of Y2 to 
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[X; Yj] is a p-submanifold, and so on. We shall denote by pn a general boundary 
defining function for a boundary hypersurface H of X. 

We now define the space M| sc . Consider in [0, e) X M X M the codimension 3 
corner C3 := {0} x dM x dM and the codimension 2 corners 

C 2 ,l ■= {0} x dM x M, C 2 ,r := {0} x M x dM, C 2fi ■= [0, 1] x dM x dM. 

We consider first the blow-up 

Ml b := [[0, 1] x M x M; C 3 , C 2 ,r, C 2 ,l, C 2 , c ] 

with blow-down map fa '■ M% b — > [0, f] x M x M. We have 7 faces on M^ b , the 
right, left, and zero faces 

rb = clos /^([0, 1] x M x dM), lb := clos ^([0, 1] x dM x M), 

zf := clos /^({O} xMx M), 
the big face bfo := /3~ 1 (C , 3) and the b-faces 

rb := clos fc 1 (C 2>i j\C 3 ), lb := clos /3 6 - 1 (C7 2 , i \C 3 ) ! bf := clos /T 1 (C 2 ,c\C 3 ). 

The closed lifted diagonal = clos 1 ( [0 , 1] x {(to, to); to e A/}) intersects the 
face bf in a p-submanifold denoted dhi^-k.b- We then define the final blow-up 

Ml sc := [M^dtfA^]. 

We obtain a set of faces denoted 2f = {zf, bf , rbo, lbo, rb, lb, bf, sc} by lifting the 
boundary hypersurfaces of Mj* and denoting sc the new face obtained by the last 
blow-up. These are shown on Figure 1 and we denote by f3 : Mj; sc — > MxMx[0, l]t 
the blow-down map, and pf a boundary defining function of f for f £ J. The lifted 
diagonal on M| sc is defined by 

A Mc := clos /? _1 ((0,1] x {(to, to), to € M}). 

Remark 2.5. The essential features of the space -M| sc are 

• for each k\ E (0, fco) the manifold M| sc PI {fc = k\} is canonically diffco- 
morphic to M S 2 C ; 

• The face zf is canonically diffeomorphic to M b ; 

• The lifted diagonal Afc iSC is a p-submanifold; 

• The vector fields d Zi lifted from the left or the right factor of M to M| sc 
are tangent to the boundary and transverse to Afc iSC except at bfo where 
they vanish simply; and 

• The face zf is disjoint from sc, reflecting the different types of asymptotics 
exhibited by the resolvent for k — and fc > 0. The intermediate face bfo 
is a transition region between the two types of asymptotics. 

Remark 2.6. This space is closely related to, but not the same as, the space in 
[4]. Either space works for asymptotically Euclidean manifolds, but the Melrose- 
Sa Barreto space used here seems to be the 'correct' one for asymptotically conic 
manifolds. 

2.2.2. Densities and half- densities. On any smooth manifold with corners X, we 
denote £lz(X) the bundle of smooth half-densities, VL b (X) the bundle of b-half 

densities. That is, SI? (X) = p~^Sl^(X) if p is a product of boundary defining 

functions for each boundary hypersurface of X. For each boundary hypersurface 

1 

H of X there is a canonical isomorphism between Q, b (X) \H and flb(H), given 
by cancelling a \dpu / Ph\^ factor at pn = 0, where pn is any boundary defining 
function for H. This does not depend on the choice of pn- Correspondingly, a 
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smooth b-half-density restricts canonically to a b-half-density on each boundary 
face H . 

We shall denote by ^ 6 1/2 (M fe 2 sc ) the lift of f\ 1/2 (M fc 2 b ) to M| sc . A smooth 

nowhere vanishing section of this bundle is v — \dgbdg' b dk / k\^ . Notice that this 
bundle is not the b-half-density bundle of M 2 SC . It can be identified with the b- 
half-density bundle except near the face sc. Indeed, as a b-half-density v vanishes to 
order n/2 at sc. We shall express our operators in terms of the bundle &l^ 2 (Mj; sc ) 

rather than Q^ 2 (M| sc ), precisely because is it generated by v (in the sense that v 
is a smooth nonvanishing section of this bundle) . 



2.2.3. Local coordinates and differential operators. Let us first consider the kernel 
of the identity operator acting on half-densities. After multiplying by the formal 
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dxdx' dydy' dk 



factor \dk/k\ x / 2 , this may be written 
(2.4) S(x - x')S(y - y') 



Let us write this in coordinates adapted to M| sc near zf, bfo and sc. The coordi- 
nates above are valid near zf provided we stay away from bfo- Near zf n bfo and in 
a neighbourhood of Afc jSC we can use x, a = x'/x,y,y' and k' = k/x'. In terms of 
these we can write the identity as 



6(a-l)S(y-y') 



dxdadydy'dn' 



XK 



which is a smooth b-half density times a distribution conormal to the diagonal 
a = l,y = y' uniformly to the boundary. 

In a neighbourhood of the interior of bfo we can use coordinates 



(2.5) 



-, k' = — r y, y' and k. 

x x' 



In these coordinates the kernel of the identity is 

dndn'dydy'dk 



5(k — n')5{y — y') 



k 



Next consider the region near bfo H sc and in a neighbourhood of Afc. sc . The 
boundary hypersurface sc is obtained from M| b by blowing up {x/k = 0, x' jx— 1 = 
0, y = y'}. The new coordinates are 



X 



x/x' - 1 
x/k 



k(- 



1 



),Y = k 



y-y 



X 



A = x/k, y, k. 



In terms of these the kernel of the identity is 

\- n / 2 5(X)5(Y) I dXdYdXd y dk I * . 

Xk 

Notice that A is a boundary defining function for sc locally near Afe iSC . We see that 
the kernel of the identity on half-densities is p sc n ^ 2 times v times a distribution 
conormal to Ak. sc uniformly to the boundary of M% sc . 

Next consider the kernel of the operator P + k 2 . We let this operator act on 
half-densities according to the formula 

(P + k 2 )(f\dg\ 1 / 2 ) = ((P + k 2 )f)\dg\ 1 / 2 . 

Another way of viewing this prescription is that derivatives are endowed with the 
flat connection on the half-density bundle that annihilates the scattering half- 
density Ifigl 1 / 2 . However, we wish to write the kernel of P + k 2 with respect to 
v = \dg\,dg' b dk/W\-l 2 . Let us therefore use the flat connection on half-densities that 
annihilates the b-half density Idgb] 1 ^ 2 - In terms of this, the kernel of P + k 2 is (in 
coordinates (x,y)) 



(2.6) 



-n/2 



(x 2 d x ) 2 -(n- 



l)x 3 d x + x 2 A S n 
= (x P b x) ® v. 



+ W + k 2 }x 



r n/2 



Notice that the vector field x 2 d x acting on the left factor lifts to M| sc to be trans- 
verse to Afe :SC except at bf where it vanishes to first order. Equivalently, the 
b-vector field xd x is transverse to Afc jSC except at sc where it has growth of order 
1 . It follows that the kernel of P + k 2 is p 2 {g times v times a distribution which is 
conormal to Ak, sc uniformly to the boundary. 
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2.2.4. Polyhomogeneous functions and index sets. Let X be a manifold with corners 
and 2f its set of boundary hypersurfaces. Let us recall that the index family consists 
of a subset £f of C x N for each f in the set 2f of boundary hypersurfaces of M% sc , 
satisfying some conditions given below. We say that a function / on a manifold with 
corners X is polyhomogeneous conormal (polyhomogeneous or phg for short) at the 
boundary with index family £ if it is conormal (i.e. if it remains in a fixed weighted 
L 2 space under repeated application of vector fields tangent to the boundary of X) 
and if for each sflwe have 

(n n (v H -z))f=o(([[ PH y) 

fffJ (z,p)e£„ s.t. Rcz<s H 

where Vh is a smooth vector field on X that takes the form Vh = Ph9 Ph + 0(p 2 H ) 
near H. This implies that / has an asymptotic expansion in powers and logarithms 
near each boundary hypersurface. In particular, near the interior of H, we have 

/= E a {z , p) p H (\ogp H T + 0(p H ) 

Z,p£&H S.t. RC2<S 

for every set, where a,( z ,p) is smooth in the interior of H (in fact, a( z ,p) will itself 
be polyhomogeneous on H). 

To make sense of this definition, index sets £f are required to have the properties 
that for each Mel the number of points (0, j) G £f with Re < M is finite, that 
if (/3, j) G £ f then ((3 + G £ f and that if j > then also (0,j - 1) G £ f . These 
conditions ensure that the space of phg functions is well-defined (independent of 
the choice of Vh) and closed under multiplication by C°°(X). We remark that the 
index family consisting of the index set N x {0} at each boundary hypersurface 
corresponds precisely to smooth functions on M 2 sc . 

Recall the operations of addition and extended union of two index sets E\ and 
E 2 , denoted E\ + E 2 and ExUE 2 respectively: 
(2.7) 

E x + E 2 = {(/?! + 02, ji +.72) I (0i, jx) G E x and (02, h) G E 2 } 
ExUE 2 = ExUE 2 U {(0,j) \ 3(0, jx) G E u (0,j 2 ) G E 2 with j = jx+32 + 1}- 
For what follows, we write q for the index set 

(2.8) {(q + n,0) I n = 0,l,2,...} 

for any q G R, and we define N and N2 to be the index sets 

(2.9) K := {(k, G No x N ; I < k}, K 2 := {(k, I)eN x N ; I < (k + l) 2 /4}- 

For any index set E and q G M, we write E > q if Re > q for all (0, j) G E and if 
(0, j) G E and Re = q implies j = 0. We write E > q if there exists e > so that 
E > q + e. Finally we say that E is integral if (0,j) G E implies that G Z. 

2.3. Operators defined on M 2 sc . We construct a parametrix G(k) for R(k) 

within a space of operators *™' £ (M; &\ /2 ) given by half-density kernels on Mj; sc . 
To each such kernels we can canonically associate an operator depending paramet- 
rically on k, as follows: we write the kernel as K eg) v where if is a distribution 
and v — Idgbdg'fjdk/k] 1 / 2 ; then (K\{ k=kl }) <S> Idgbdg'^ 1 / 2 defines an operator on half- 
densities on M for every kx G (0, fco). For example, the resolvent of the Laplacian 
on R 3 is represented by 



1 e ~ fc l z - z 'l , / lf 2 ,dk,x/2 {xx')- 3 / 2 e - k \ z - z 



v 



(2.10) ~ ®\dzdz'\ 1 / 2 ®\^ ] 

At: \z - z'\ 1 1 k 1 47r \z - z'\ 

We want a 'variable coefficient' class of operators modelled on the free resolvents 
on R™. The essential properties we shall require of our kernels are (i) conormality 
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at the diagonal, uniformly up to the boundary, and (ii) polyhomogeneity at the 
boundary away from the diagonal. We remark that the diagonal Afc iSC C M 2 sc is 
a p-submanifold transverse to each boundary hypersurface so conormality makes 
perfect sense here. 

The right-hand expression in (|2.10p for the resolvent on R 3 shows that the kernel 
vanishes to order —2 at bfo and order — n/2 at sc. Moreover, the kernel vanishes to 
infinite order at bf, lb and rb. This motivates the following definitions: 

Definition 2.7. The class ^™' £ (M; £l^ 2 ) of operators on half-densities on M is the 
set of half-density kernels K on M 2 sc that can be written as a sum of two terms 
K = K % + K 2 , where 

(i) psl" 2 K\ is supported near and conormal to Afc jSC of order m, uniformly up to 
the boundary, and 

(ii) Psl 2 K% is smooth in the interior of Af| sc and classical conormal at the 
boundary with respect to some index family £, i.e. an index set £f at each boundary 
hypersurface f. Moreover, we assume that 

(2.11) £bf, £ r b arid £ib are empty 
and that (in the notation of (|2.8[> ) 

(2.12) £ z f, £bf and £ sc contain the C°° index set 0. 

Definition 2.8. The class \j™'( abf o' azf ' asc ) ,e (_^f. Q^ 2 ) of operators on half-densities 
on M is the class of half-density kernels K that can be written as a sum of two terms 
K = K 1 + K 2 , where K 1 = pl¥ ° p"f P%°° K i and K i and K ? are in Definition 
with (12.12[) replaced by 

(2.13) £bf 07 £zf and £ sc contain the index set abf , a z f and a sc respectively. 

Remark 2.9. The order —n/2 of growth of the kernel at sc reflects the fact that we 
are expressing the kernel in terms of b-half-densities when scattering half-densities 
would be more natural here. 

Proposition 2.10. Suppose that 

A E ^ a ^ a ^ A (M -hl /2 ) and B G *f ' (6 **' 6 * 6 - )iB (M; fij /2 ) 
satisfy (|2.1ip and (|2.13p . Then the composition A o B is an element of 

^m+m' ,{abf +bbf ,a zS +b zS ,a sc +b sc );<Z fj 1 / 2 ) 

with 

Gsc "Asc ~t~ sc 

&zf = (A z f+ ¥i z f)u[A r b + "Bib ) 

e &/o = { A ib + ' B rb )n(yib /o + 3 6 / ) 

£ib = {Am + T> zJ )u(A bfo + ■BzbJ 
Grb = (A zf + T, rba )u(A rba + B 6/o ) 

Proof. See Section H3 □ 

Corollary 2.11. Suppose E G y™^ a Vo> a *f> a ^ E ( M . q 1 / 2 ) w h ere m < 0, each of 
<^bf Q ,o JZ f,a sc is positive, and £/ > for all f. Then for large enough N, E N is 
Hilbert- Schmidt for each k > with Hilbert- Schmidt norm \\E N {k)\\Hs ~ * as 
k — * 0. In particular, E is compact for N large enough, uniformly (in the sense 
above) as k — > 0. 



(2.14) 
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Remark 2.12. The operator E will arise as the error term in a parametrix construc- 
tion, i.e. we will have 

(P + k 2 )G(k) = ld+E(k) 
where G(k) and E(k) are kernels on M% sc . We shall say that E = E(k) 'iter- 
ates away' if it has the properties in the corollary above. The point is that in 
the Neumann series for (Id +E) -1 , each term becomes sucessively milder, and the 
Neumann series can be summed asymptotically modulo a smooth kernel on M| sc 
that vanishes to infinite order at the boundary. Moreover a finite number of terms 
in the Neumann series gives an error that is compact uniformly in k. 

Proof. Choose S > so that —m > 5, dbf , a zf> &sc > 8, and £f > S for all f. Then, 
by Proposition 12.101 E N has a similar property with respect to NS. For NS > n 
the kernel is continuous and vanishes to order greater than n at each boundary 
hypersurface of M 2 sc which implies that it is Hilbert-Schmidt for each k > with 
Hilbert-Schmidt norm tending to zero as k — > 0. □ 

2.4. Normal operators. Assume that 

(2.15) £ zf , £ b f and £ sc > 0. 

Then the restriction of A G 4 rm ' £ (M) to each of these faces is well-defined and 
can be viewed as an operator on functions on a model space. On zf this model 
operator is a b-pseudodifferential operatoid acting on half-densities on M, on bf 
it is a pseudodifferential operator acting on half-densities defined on S" 1-1 x (0, oo) 
and on sc it is a family of convolution pseudodifferential operators, parametrized 
by z G S*" -1 and by k G (0, fen], acting on half-densities defined on R". We call 
the model operator at f G {zf, bfo, sc} the normal operator of A at f and denote it 
HA). 

Proposition 2.13. The normal operators respect composition: if A G x f , " l " A (M, ) 
and B G # m ' S (M, f^ /2 ) satisfy (j2~T5)) and 

(2.16) A rbo + S Z6o > and A Wa + 'B rbo > 
then 

I zf (AoB)=I zf (A)oI zf (B) 

(2.17) I bfo (AoB)=I bfo (A)oI bfo (B) 

Isc ( A oB) = I SC (A) o I SC (B) 

Proof. The third composition property is just that of the scattering calculus with 
a smooth parameter k. 

To prove the first identity, note that I z f(A o B) is the pushforward of ir* L A (g> 
TtftB ® -k* c v ® \dklk\~ 1 ! 2 restricted to the boundary hypersurface of M^ sc which 
is the lift of M 3 x {0}. This boundary hypersurface is canonically isomorphic to 
M b and the map to zf is the canonical b-fibration. This implements the product 
in the b-calculus, so the first line follows immediately. (We remark that there is 
another boundary hypersurface of sc that projects under ire to zf but under the 
assumption (|2.16|) it does not contribute to the leading order term.) 

The second identity follows similarly. □ 

Remark 2.14. Following on from the last point in Remark 12.51 the face bfo is the 
space (dMxI) 2 , where I is the compactification of the interval (0, oo) K parametrized 
by k = k/x, with a 'b'-blowup at k = k! = and a scattering blowup at n = k' = oo. 

■^This is not quite true, since there may be a poly homogeneous expansion of I z i(A) at the 
boundary of zf which is not allowed in the b-calculus as defined in However this is an 

inessential point. 
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Looking ahead to Section [31 we see from (|3.2p and (|3.8[) that the normal operator 
of P at bfo is of b-type at k — and of scattering type at K = oo, so we see the 
transition between b- and scattering type behaviour quite explicitly. 

2.5. Strategy. We attempt to solve the equation 

(2.18) (P + k 2 )G{k) = Kw 

where G — G(k) is in P^ ^~ 2 ' E (M 2 sc ;Cl^ 2 ). Of course we cannot do this in one 
step, so we first try to write down a parametrix G with an error term E = E{k) 
which iterates away in the sense of Remark 12.121 

In the simplest case that P has no I? null space nor a zero-resonance, our 
parametrix G will have the index set at zf and the index set —2 at bfo- Let us 
write G° f for the restriction of G to zf (which is canonically defined since G is a 
b-half-density) and G^ 2 for the restriction of k 2 G to bfo. In general, if we need 
to consider more terms in the Taylor series of G we shall write G J f , j £ N, for the 
jth term in the Taylor series of G at any face f. Of course we need to specify a 
boundary defining function for f if j — 0. We shall always use the boundary defining 
function k at zf rbo, bf or Ibo; this is the most convenient choice since it commutes 
with the operator P + k 2 . We remark that k is only a boundary defining function 
in the interior of each of these faces, but this is of no importance; we only have to 
remember that the coefficients Gj then may become more and more singular at the 
boundary of f as j increases. 

With this notation, then we see that if G £ p^ 2 ^^ 2 ^ (Mj 2 sc ; &l^ 2 ) then our 
requirements are that E® { , E^ { and E® c all vanish. This amounts to inverting 
model elliptic operators at each of these three boundary hypersurfaces. It turns 
out that in high dimensions, n > 5, this is sufficient to obtain an error term that 
iterates away. In low dimensions or in the presence of L 2 null space or a zero- 
resonance we need to work harder, that is, to solve further model problems. Wc 
shall begin with the simplest case first. 

3. Resolvent kernel when (M,g) is asymptotically Euclidean 

Let M, g) be an asymptotically Euclidean manifold. We shall construct a para- 
metrix for (P + k 2 )^ 1 with an error that 'iterates away' in the sense of Remark l2.12l 
Throughout this paper we assume that P has neither nontrivial L 2 null space, nor a 
zero-resonance. Initially we shall assume that n > 5, as the parametrix construction 
then is as straightforward as possible. In section [3771 we treat the cases n = 3, 4. 

3.1. Singularity at the diagonal A& iSC . The operator P + k 2 is elliptic as an 
element of v[< 2 >( 2 > ' );*(M, f^ 2 ), i.e. its symbol multiplied by p^ 2 Q is elliptic uni- 
formly on iV*(Afe jSC ). We specify the symbol of G(k) to be the inverse (in the sense 
of the noncommutative product corresponding to operator composition) of that of 
P + k 2 modulo symbols of order — oo. This is consistent with membership of G{k) 

in v]/- 2 ;(- 2 > > )'*(M, fi^ 2 ) and determines the symbol of G(k) modulo symbols of 
order — oo. 

3.2. Term at sc. The normal operator of P + k 2 at sc is the family of operators 
A + k 2 acting on half-densities on R™, parametrized by z £ S*™ -1 and k £ (0, fco]. 
We specify the normal operator I sc (G(k)) to be the inverse operator (A + fc 2 ) -1 
which is well-defined for k > 0. 
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3.3. Term at zf. Let us start with the face zf, which is canonically identified with 
the b-blow-up M| := [M x M; dM x dM]. On this face zf, the equation is really a 
b-elliptic equation in disguise: we see from (|2.6|) that P = xP\,x (where we interpret 
P and Pb as acting on half-densities, and P annihilates the scattering half-density 
while Pfc annihilates the b-half-density) where the operator p, is self-adjoint with 
respect to the b-metric <?& and elliptic in the b-calculus. Moreover, its normal 
operator iff (p,) is invertible. Therefore, by the theory of b-elliptic operators given 
by Melrose [HI Sec. 5.26], there is a generalized inverse Qb, a b-pseudodifferential 
operator of order —2 for the operator Pb on L 2 , such that 

(3.1) p b Qb = Q b p b = id -n 6 . 

where Tib is orthogonal projection on the L 2 kernel of Pb- From .Hi], an element tp 

i 

in the range of lib is in x^^ 1 L°° (M, ) (see Theorem 12.11 and Proposition 12. 3|) . 

for all e > 0. This implies that x~ x tp G x n - 2 L 2 (M, fijjL) satisfies P(x"V) = 0. So 
x~ X Lp is either a zero-mode or a zero- resonance. By our assumptions on P, we have 

n 6 = 0. 

Again from [TTj . the kernel Qb is conormal at the b-diagonal Afc !SC nzf, uniformly up 
to the front face zfnbfo, and is the sum of three terms: one is supported close to and 
conormal at the b-diagonal, uniformly to the boundary; one is polyhomogeneous on 
M 2 with the index set £(Qb) = (£ff (Qb), £ r b(Q&), £ib(Qb)) where rb,lb,ff denote 
the right boundary, left boundary and front face of the blow-up M 2 , and 

TL 

£rb(Q 6 ) - £ib(Q 6 ) = {(n/2 - 1 + I, k); < k < I e No} = - - 1 + N, £ ff = 

and the other term is the lift of a polyhomogeneous distribution on M x M, with 
index sets £ r b(Qh) at the right and left boundaries. 
So we have 

PbQb = Id . 

Therefore we set 

c° f = (xx'r'Qb. 

3.4. Leading term at bfo We use coordinates (|2.5p near the interior of bfo in 
terms of which we have int(bfo) = {k = 0, k, k' G (0, oo)}. The operator P + k 2 in 
these coordinates reads 

P + k 2 = /c 2 K"t-^ _ (nd K ) 2 + Agn-i + (n- 2) 2 /4 + k 2 +W^k%~ 1 . 

Since VF G xDiffb(M), the first germ of P + k 2 at bfo is at order 2, we consider 

hf (k- 2 (P + k 2 )) = ( - (nd K ) 2 + A s »-i + (n - 2) 2 /4 + k 2 )^- 1 . 

Here we wrote P with respect to the flat connection on half-densities annihilating 
Icigl 1 / 2 . In terms of the b-flat connection annihilating the b-half-density {dgbl 1 ^ 2 , 
our operator is 

(3.2) ^ — (nd K ) 2 + Ag«-i + {n-2) 2 /A + k 2 ^- 1 := k^Pm^" 1 . 

Setting Mo := (0, oo) x S n , we can view int(bfo) as Mq x Mo- The operator Pbf 
acting on the left on Mo x Mq on b half-densities /(k, y)|K -1 dK<i?/| 2 has an inverse 
in terms of spherical harmonics {<frj{y))j from Sturm-Liouville theory: if H be the 
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Heaviside function, I W) K U the modified Bessel functions, then 

oo 

g bfo :=E n ^ (/i+t-i(«)^i-i+i(« , )-ff(«' - «) 



(3.3) 



dndydn' dy' 



solves on b-half densities on Mq x Mo 

-Pbf Qbfo = Id. 

Therefore we set 



(3.4) 



G 



bin 



(kk')Q 



bin 



3.5. Consistency. Before proceeding further we need to check that our specified 
values for a(G(k)), G® c , G° f and G^f are consistent. 

We first check consistency between G" f and G^[ . To do this we must verify that 
the restrictions of p 2 ft) G° f and p 2 f G 2 fo agree at zf n bf , where p z f ■ pbf = k. It 
is most convenient to take pj^f = {xx')~ x and p 2 f = kk'. Thus, we compare the 
restrictions of (xx')G^f and (kk') _1 G|^ at the intersection bfoflzf. From G° f we get 
the normal operator (in the sense of the b-calculus) of Qb- This is the inverse of the 
normal operator of which has the exact expression (sD s ) 2 + (n/2 — l) 2 + Ag»-i 
where s := (3* (x/x') S K is a coordinate on the fibers of ff. The solution is given by 
standard Sturm-Liouville theory in terms of the eigendecomposition of Ajn-i by 

~ e -C7+f-l)|log*| 

(3.5) E n ^ ■>.•■.. o \dydy'da/8\*. 



2 j + n - 2 



On the other hand, to compute the restriction of (Kn')G h 2 o at the intersection 
bfo H zf, we recall the definitions of I v , K v and their asymptotics 

7T J„(-z) 
2 si 



(3.6) I v {z) 



fc!i> 



i 



z_ 



while for = rt G Z, /„ is entire, /„ = i_ n and -ftT,, is given by the limiting value of 
the right hand side as v — > n. As z — » 0, we have 

1 



(3.7) 



i> 



r(i/) ,z 

2 V ' 



j/ > 0. 



Also K v {z) = 0(e~ c ' 2 ') for some c > as \z\ — > oo with z,i/£l. At bfoflzf we have 
k = k' = and only the leading behaviour (|3.7p of I u and -K^ at k — contributes 
to this restriction. We see from this that the restrictions agree identically. 

We next check consistency between G® c and G^ 2 g . The face sc is fibred over 
(fc, y) with fibres having a Euclidean structure. Let z be a Euclidean variable on 
the fibres of the scattering double space M 2 C ; then kz is a Euclidean variable on 
the fibres of sc. We have defined G® c — k n ~ 2 f(kz)\dz\, where f(z) is the kernel of 
(A + I)" 1 on R n and hence k n ~ 2 f(kz) is the kernel of (A + k 2 )- 1 on 1". 

At bfo, let us write the normal operator G^ 2 in terms of the scattering half- 
density k n \dgdg' dk/k^/ 2 which restricts to the half-density \n n ~ 1 dndyK' n ~ 1 dK'dy / \ 1 / 2 
at bfo- That is, let 

^Mo = Qbi k n \dgdg' dk / k] 1 / 2 . 
Correspondingly we use the flat connection on half-densities that annihilates the 
scattering half-density. The PDE for Qbf then takes the form 



(3.8) 



2 Ac 



K K 
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This is just the equation A + 1 on Mq viewed as R n in polar coordinates (k, y), with 
k playing the role of the radial variable. It follows that we have Qm — f{ K y~ K 'y')i 
where we view ny and n'y' as lying in W 1 . 

The boundary hypersurface bfo may be obtained from Mq by first compactifiying 
Mo, then forming the b-stretched product, and then blowing up the boundary of the 
b-diagonal at the 'oo' end only. This last blowup is at scflbfo and this face is fibred 
over S n ~ l with fibres that are Euclidean with Euclidean variable z — ny — n'y' . It 
follows that the restriction of Qbf to this face is f(z). We finally observe that z 
agrees with the limiting value of the coordinate kz taken from sc. This shows the 
compatibility between G® c and G^f 2 • 

Next we check consistency between the symbol of G(k) and the normal operators. 
This is immediate because each of the normal operators at sc, zf and bfo solve an 
elliptic pseudodifferential equation so the symbol at the diagonal on each face is 
determined uniquely modulo symbols of order — oo, and hence necessarily agree 
modulo symbols of order — oo with that coming from the interior of Afe jSC . 

We next consider the behaviour of these normal operators at the remaining 
boundary faces. We note that G® c decays to infinity order at the boundary sc n bf, 
uniformly down to k = 0. (This is because a boundary defining function for bf near 
sc n bf is k\z — z'\.) Similarly, the exponential decrease of K v (k) as K — » oo shows 
that Gy* decays rapidly at the intersection bfo fl bf, bfo n lb and bfo Pi rb. Finally, 
both G^ 2 and G° f are polyhomogeneous at lbo and rbo with index set n/2 — 2 + N. 
Moreover G^ 2 q has index set at zf given by 

(3.9) if n is odd, + {(n - 2 + 21, 1) | / e N } if n is even. 

We conclude that there exists a G = G(k) 6 *- 2 >(- 2 >°>°); s (M, fi^ /2 ) satisfying 
9 SC = 0, Sbt = -2, Sib = Srb = n/2 - 2 + X, S zf is given by the index set 

at rb, lb,bf is empty, and G has G^ 2 , G^ f , G°sc as leading term at bfo,zf, sc, with 
symbol at the diagonal given in 13.11 

3.6. Exact resolvent. Applying P + k 2 to our parametrix G(fc), we obtain 

(P + k 2 )G(k) = ld+E(k) 

where E(k) e $-°°,(M,i);£( M; qV^ where g sa ti s fies £ zf > 1, £ b f > 1, £ sc > 1, 
£ib„ > n /2, £ r b > n/2 — 2 and (|2.1ip . (We emphasize that since the operator 
P + k 2 vanishes to order 2 at bfo, we automatically gain two orders at bfo, and 
the fact that we solved the model problem at bfo gives us the additional order of 
vanishing. Hence E(k) is three orders better at bfo than G(k).) In particular, 
£f > for each boundary hypersurface /. It follows from Corollary 12.111 that E(k) 
is Hilbert-Schmidt with Hilbert-Schmidt norm going to zero as k — > 0. It follows 
that ld+E(k) is invertible for small k. Moreover, E(k) 'iterates away' in the sense 
of Remark 12.121 It follows from standard arguments that (ld+E(k))~ 1 lies in the 
calculus. Inspection of (j2. 14(1 and an inductive argument shows that for I > 2 

E(k) 1 =e *-~' (w);e! (M:^ /2 ) 

where £ ( is the index set 

4 = i, £i f ,£ bfo cl + N|, 

E l tb0 C~ + l-3 + 7ii, fi^c^ + i-l + N, 
with Ki := {(k,j) e N x N ; j < l(k + 1)}. Let us write 

(Id+Eik))- 1 = ld+S(k). 
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Then S(k) is an element of (M, h\ /2 ) for some index set 8 determined 

by £. Inspection of (|2.14p and an inductive argument shows that 

S zf ,Sbf C 1U(2 + K 2 ), S lbo c(|+^ 2 ), S rbo C (|-2)U(|-1 + N 2 ), S sc = l 
where N 2 is defined by 

N 2 := {(k,j) e No x N ; j < (fc + 2) 2 /4}. 

It follows that the resolvent kernel can be written R(k) = G(k) + G(k)S(k) and 
also lies in the calculus. A final application of the pushforward theorem gives (1 1 . 3[) . 
with 

Krf cOU(l + N 2 ), ^ib ,3?rb C (- -2 + N 2 ), 
(3.10) * 

31bf C -2U(-1 + K 2 ), ^ sc = 0. 

Moreover the correction term G(k)S(k) vanishes to first order at zf and bfo so 
our parametrix is exact to leading order at these faces. Moreover, the asymptotic 
expansion of the resolvent R(k) as k — > can, in principle, be generated to any 
given order at each boundary hypersurface from computing a finite number of terms 

in the Neumann series, i.e. from G(k){ld-E(k) + E(k) 2 h (-l) N E(k) N ) for 

sufficiently large TV. 

Remark 3.1. We emphasize that the index families given above are not sharp as 
far as the logarithmic exponents are concerned. 

3.7. Leading term at rt>o and lbo. We can improve the accuracy of our para- 
metrix by specifying G™^ 2 2 and G^ 2 2 appropriately. When this is done correctly 
we obtain a better error term E{k) with leading behaviour at rbo and lbo at order 
n/2 — 1 instead of n/2 — 2. This is useful for two reasons. The first is that this 
allows us to construct the resolvent in dimensions 3 and 4. The second is that it 
allows us to determine the leading behaviour of the resolvent itself at rbo and lbo, 
which turns out to be crucial in understanding the L p boundedness properties of 
the Riesz transform. 

The term G™/ 2 2 must solve the equation P^xG™/ 2 2 ) = and match correctly 
with G° f and G^ 2 . Let us first consider the behaviour of G° f = [xx')~ x Qb near 
zf n rbo ■ Localizing near rbo > we have by Theorem [ 



(3.11) Q b = (v(z)x'^ 1 + 0(x' S )) 



dxdx' dydy' 



for some v{z) with a polyhomogeneous expansion in z — (x,y) as x — > 0. In order 
to match with the normal operator at the front face (here zf n bfo) v can grow no 
faster than x ~ n / 2+1 as x — > 0. Proposition ^. 31 shows that there is a unique element 
of the null space of P& (up to scaling) with this property, which has the behaviour 
cx -n/2+i _j_ q^ x — n/2+2^ ag x q wnere c j s constant in y. Matching with the 
normal operator of Qt then gives the leading asymptotic 

1 

(n - 2)Vol(S'"- ] ' 
To match with G° f then we take 



"CO = 7 KyTt uc n -i\ x ~* +1 + 0(x~%+ 2 ). 



G^- 2 =x-i V (z)f( K >,y>) 



dxdx' dydy' 



where /(«/, y) ~ (r(n/2 — l)2 n / 2 2 ) 1 k' n ^ 2 2 as k' —> and otherwise is deter 

'b{ 



mined by the matching condition with G bf 2 . At the intersection rbo H bfo, n = 



18 



COLIN GUILLARMOU AND ANDREW HASSELL 



and k' £ (0,oo). In particular we have K < k' . Examining (|3.3p and (|3.4p we see 
that the leading behaviour of G^ 2 at bfo n rbo is 



1 



T(n/2 - 1) 



2 K K n /2-l(K ). 



Hence we set f{n',y) = (T(n/2 - l))- 1 2- n ^ +2 K'K n/2 _ 1 (K r ), i.e. 



g: 



1/2-2 _ x 1 v{z)K'K n /2-i(n') 
b ° ~ r(n/2 - l)2»/2-2 



and this matches with both G° f and G bf 2 . 



We now set 



dxdx'dydy' 



xx' 



G 



i/2-2 _ ^ «(g0^gn/2-lW 



lb„ 



T(n/2 - 1)2™/ 2 - 



dxdx'dydy' 



xx' 



to satisfy formal self-adjointness at this order. We find that the error term now is 



,(i,i,i);£'( Mj ^V 2 ) ; where £ / sat i snes f or some £^> 1- 



e, £ bfo >l' 



#(fc) e 

£s C > 1 - e, £; bo > n/2 + 1 - e, £J. bn > n/2 - 1 - e and (j2~TTj) for some small 
e < 1. Then one can check that G(k)(E(k)) N has an index set at rbo which is 
always > n/2 — 1 — e. That means that the resolvent kernel itself has an index 
set at rbo,lbo included in (n/2 — 2) U Z for some index set Z > n/2 — 1 — e, and 
the term of order n/2 — 2 satisfies i2^ 2 ~ 2 = G™/ 2-2 , i.e. we have found the ac- 
tual leading behaviour of the resolvent kernel at rbo- Moreover, the construction 
now goes through in dimensions 3 and 4, assuming that the L 2 null space is triv- 
ial and that there are no zero-resonances, since now the error term E(k) vanishes 
to a positive order n/2 — 1 at the right boundary as opposed to n/2 — 2 as previously. 



Remark 3.2. We recall that v(z) solves Px% ~ 1 v = and x^^ 1 v(z) is asymptotically 
constant at infinity. If M has only one end, then Theorem 12.11 implies that this is 
the only function (up to scaling) satisfying these properties. Hence, if M has only 
one end and V = 0, then x^~ 1 v(z) is necessarily constant. 



4. Resolvent kernel when (M,g) is asymptotically conic 

We essentially give the differences with the asymptotically Euclidean case. First 
recall N a := {(n - 2)/2 = v < v x < . . . } is defined in ([2~3]) . 

4.1. Singularity at the diagonal Afe iSC . This is completely similar to the asymp- 
totically Euclidean case. 

4.2. Term at sc. The fibers of sc are Euclidean spaces radially compactified and 
the normal operator of P + k 2 at sc is the family of flat Laplacians A + k 2 acting on 
half-densities on these Euclidean spaces, parametrized by y' € dM and k € (0, fco]. 
Indeed if y%, . . . ,y n are local coordinates on dM, the functions zq — k/x—k/x', Zi := 
k(yi/x — y'Jx') give Euclidean coordinates on the fibers sc and the normal operator 
of A g + k 2 at the fiber of sc with basis point (fe, y') at sc is A: 2 (A 2 + 1) where A z 
is the Euclidean Laplacian — d 2 — ^2n {y')d Zi d z ., ho(y') being the Euclidean 
metric on this fiber induced by fiQ. We specify the normal operator I sc (G(k)) to be, 
in each fiber, the inverse operator k~ 2 (A z + 1) _1 which is well-defined for k > 0. 
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4.3. Term at zf. As in the case dM — S"™" 1 , the theory of b-elliptic operators 
given by Melrose [HI Sec. 5.26] shows that there is a generalized inverse Qb, a 
b-pseudodifferential operator of order —2 for the operator Pb on L 2 , such that 

P b Qb = QbPb = w -n 6 . 

where Ilf, is orthogonal projection on the L\ kernel of Pf,. As in Section [3 due to 
our assumptions on P, we have 

n h = 0. 

As in Sect ion [331 Qb is conormal at the b-diagonal Afc jSC nzf, uniformly up to zfnbf , 
and is the sum of three terms as described there, except that the polyhomogeneous 
term has the index set £(<9&) = (£ff(Q&), £rt>(<2&), £it>(<2i>)) where 

£rb(Q 6 ) - £ib(Q 6 ) c {(u, k) e R x N ; v e N, k + l < e N o; n < v}}, 

£ ff = N x 0. 

So we have PbQb = Id and we set G° f = (xx')~ 1 Qb- 

4.4. Leading term at bfo. This is essentially similar to the case dX — S 71 ^ 1 
except that we have to replace Agn-i by Aqm and Pbf on M — (0, oo) x dM has 
an inverse given in terms of eigendecomposition of Aqm'- if E Vj :— ker(AgM ~ v] + 
(n - 2) 2 /4) where ij G N 9 , that is 

OO 

Q bfo := U Ei (iu, {*)K Vj {k')H{k' -k) + I Vj {n')K Vj (k)H(k - k')) 

dudydu'dy' 2 
kk' 

solves -Pbf„Qbf = Id on b-half densities on Mo x Mo. The index set of Qm at zf is 
of the form + Z where Z is an index set satisfying Z > min(2, n — 2) — e for any 
e > 0. Therefore we set 

G bf„ = (*«0Qbf o < 

Remark 4.1. If M° has j >2 ends, i.e. <9M has j > 2 components S±,...,Sj, then 
bfo has j 2 components. In this case, all of the off-diagonal components of G^ f 2 are 
zero, while the diagonal components are completely independent: the ith diagonal 
component only depends on the metric h(0) at Si. 

4.5. Leading term at rbo and lbo. The leading terms G" b /2 - 2 and G™ /2 " : are 
defined similarly to the previous section since Vq = (n— 2)/2 here as in the asymp- 
totically Euclidean case. However, we remark on what happens when there is more 
than one end, say j > 2 ends. Then from matching with G^ f 2 we obtain bound- 
ary conditions for the function v(z) from Section f3. 71 on the ith component of rbo: 
namely, it should tend to a nonzero limit at Si and to zero at all other ends (cf. 
Remark 14. ip . Again Theorem 12.11 implies that there is a unique function (up to 
scaling) with this property. 

This is sufficient to deal with dimension n > 4 since we gain an error that has a 
positive order of vanishing at rbo, lbo. However, to deal with dimension n — 3, it 
is not sufficient and we have to construct the terms G 1 ? 1 such that < za- < 1/2 

rbo J — ' 

in the parametrix. Actually it is possible to compute the actual second asymptotic 
term at rbo m any dimension by this method, when V\ — vq < 1, so we will write it 
for completeness and because it will imply an interesting result for Ricsz transform. 

For that, we notice that Qb has asymptotic expansion near rbo 
Q b = £ x^v Vj {z,y') + O{x' V0+1+e ) 

Vo<Vj<V + l 
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for some e > and PbV Uj (z, y') = with v Vj growing like x Vj . Checking consistency 
at bfo with Is(Qb) shows that 

(4.2) v Vj {x,y,y') = x -"i + 0{x~ v ^ logs). 

Now it suffices to set 



which matches with Gj^ and G° f . 

4.6. Consistency. This is almost the same as in Section [3 we just have to change 
the normal operator of Qb on zf n bfo since the normal operator of Pb has the 
expression ©jeN (( S D S ) 2 + Vj) (where s :— (5*{x/x') 6 1 is a coordinate on the 
fibers of ff) after decomposing according to the eigendecomposition of Aqm- The 
solution is given by standard Sturm-Liouville theory 

~ e -«i|logs| A 

2^ n Bj . — — \dydy'ds/s\ = . 

j'=o Uj 

Then everything is similar to the case dM = S 71 ^ 1 except that G^"/ and G° f 
are polyhomogeneous at lbo and rbo with another index set £ r b(Q&) — I > n/2 — 2 
which does not necessarily contain half integers. 

Therefore, there exists a G = G(fc) 6 ^r-2,(-2,o,o) ; S( M) qV 2 ) satisfying S sc = 0, 
Szf C U Z, 9 b f = -2, Sib = Srb C U Uo <„ j<Vo+1 (vj,0) Uvq + Z for some index 
set Z > — e for any e > 0, and G(k) with the values at each face prescribed by our 
models. 

4.7. Exact resolvent. Applying P + k 2 to G(k), we define the error E(k) by 

(P + k 2 )G{k) = ld+E(k) 

where E(k) E ^-^'^'^(M,^ 2 ), and £ satisfies £ zf > 1 - e for any e > 0, 
£bf > 1) £sc > 1, £ib > n /2, £rb > i/2 — 1 — c for any e > 0. In particular, 
£f > for each boundary hypersurface / and Corollary 12 . 1 1 1 shows that the method 
detailed for dM — S 71 ^ 1 applies, and the resolvent R(k) has the same properties 
except the index set which becomes 

Xi C (0, 0) U (1 + T), % ho , % h(> C |J (vj - 1) U (n/2 -l + T), 

(4.3) v <u } <u +l 

tt b f C (-2,Q)u(-l + 3-), ft sc = 

for some index set £F > — e for any e > 0. Moreover Proposition 12.101 shows that 
the first asymptotic terms i?^ b 1 of the resolvent at rbo are given by G^ b 1 for 
Vj < vq + 1 = n/2 — 1. 

Remark 4.2. Suppose that (M°,g) is a metric cone. Let Ao be the Laplacian on M° , 
defined as the Friedrichs extension of the operator with domain given by smooth 
functions with support compact and disjoint from the cone point. (This operator 
is essentially self-adjoint in dimensions n > 4 but not for n = 2, 3.) The resolvent 
i?o = (Ao + k 2 )~ 1 scales as A~ 2 under the transformation x — » Ax, x' — > Xx' 
and k — ► Xk and is therefore determined by its leading term at bfo, (Ro)^f - It is 
not hard to see, by expanding in eigenfunctions on the cross-section dM, that this 
model is given precisely by (|4.1[) . Now choose a cutoff function <j> that is equal to 1 
near infinity and near the cone point. Then passing to 4>Rq4> removes singularities 
of the kernel of Rq due to the cone point, and we see that this kernel satisfies all 
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the properties of the resolvent listed above. We shall use this fact in the following 
section. 

5. RlESZ TRANSFORM 

5.1. Boundedness on L p of a class of operators. We first determine a class of 
operators in our calculus which are bounded on L P (M). 

Proposition 5.1. Let a, /3, <5 > 0. Assume that A = A(k) G ^;}~ lfi ' 0) ' A (M, Q* /2 ) 
with index family A satisfying 



Aw = A rb = A bf = 0, A bfo > -1, A sc > 0, A zf > -1, 
Then J*° A{k) dk is bounded on L P (M) fo 



Ti 7~l 

Arb >^-l-a, A ibo >--l-[3. 



n 

< p < 



n- (3 

Proof. Let us first set 8 > so that A z f > S and Abt > S. We give the details 
when the kernel is localized to the region x < e,x' < e. When the kernel of A is 
localized to the region x > e/2, x' > e/2 then A(k) is a bounded family of properly 
supported pseudodifferential operators of order —1 and the result is classical, while 
if we localize the kernel to x > e, x' < e/2 or to x 1 > e, x < e/2 then the argument 
is similar (but simpler) to that below for the kernels A\ and A^. 

In the region x < e, x' < ewe can use coordinates (x, y) in a collar neighbourhood 
of the boundary of M. We define s — x/x' and write dy for the measure density 
on dM. 

Let us decompose A = A\ + A p h g where A\ is supported close to the diagonal 
Afe jSC C Mf sc and A p h g is polyhomogeneous on M| sc with the same index family 
as A. As above we assume both kernels are supported in x < e, x' < e. We further 
decompose A p h g = A\ + A^ + A3, where A\ is supported away from sc and in 
s > 1, A2 is supported away from sc and in s < 1 and A3 is localized near sc. We 
decompose A a = A4 + A 5 where A4 is localized in k > x' and A 5 in k < x' . Note 
that we can assume that l/2<s<2on the support of A 4 and A 5 . 

It is natural here to switch back to writing the kernel of A as a multiple of the 
scattering half-density \dgdg'\ x ^ 2 , since our L P (M) spaces are with respect to the 
Riemannian measure \dg\. Indeed we ignore the half-density factors from here on 
(which are not natural for p ^ 2). Then the order of vanishing of the kernel is 
n — 1 + e at bfo for some e > 0, n - 1 - a at rbo and n — 1 — [3 at lb , and at sc 
and zf, due to the relation \dgbdg' b \ 1/2 = (xx') n l' 1 \dgdg'\^l' 1 . 

Consider the kernel A\ . Since it is supported away from sc we may use boundary 
defining functions x' / (x' + k) for rb, x + x' + k for bf and (x' + k)/(x + x' + k) for 
rbo. So we can bound the kernel of J A\dk by 

< C N x' n+5 s a+s / k- 1+s (l + ky N+n dk < C N x' n+5 s a+s 
Jo 

where N > is arbitrary and we used change of variable kx' = k. We look at L p 
boundedness < 
need to show 



boundedness of J Q A\{k)dk. Let f(x,y) be a function on [0, e] x x dM. Then we 



Ho IIo 



f(x',y')s a+s x'- 1+S dx'dy 



x-^dxdyKCWfWl 



Lp(M)- 
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We bound the left hand side using Holder's equality (with q = p/(p — 1)) to obtain 





J Jo J.C f[r '-"' lr ' : ''''''' 

£ II! ( /f |/( *' y,) ^)( ff^- 1 -*™''^*)*^^ 

< C\\f\\ p p J j\ p( - 1 ~ a)+n+1+{p - 1) x- n - 1+ ^ a+s Uxdy 

= C\\f\\ p p Jj\v 8 - l dxdy<C'\\f\\l. 

where we used p < n/a to integrate the x' integrand of the second line. 

The same argument applied to the formal adjoint of A 2 shows that A 2 is bounded 
on L p for p > n/(n — (3). 

For the kernel A3, we may localize near a fiber of sc with basis point y', then if 
2/1, . . . , y n are local coordinates near y', we can use the boundary defining function 
(l + fc 2 |z — z'| 2 ) -1 / 2 for bf, where zq = 1/x— l/x 1 and z% — yi/x — yi/x 1 ', we also use 
x + x' + k for boundary defining function of bfo. However, since A3 is by assumption 
supported near sc, we may assume that k > (x + x 1 ) on the support on A 3 and 
hence we may take k as a boundary defining function for bf . This allows us to 
bound J A 3 (k) dk by 



k ° 1 N 



i ± + k2\z 

for any N. By changing variable to k\z — z'\ in the integral we see that this is 
bounded by min(l, \z— z\~ n ~ e ), which shows that it is bounded on LP for 1 < p < 00. 

For the kernel A4, we use same notations for z and use coordinates W := k(z — 
z'), y, = x'/k, k as coordinates. Then W = defines the diagonal. Since A4 is 
conormal to the diagonal of conormal order — 1, and vanishes to order n — 1 + S at 
bfo, the kernel is bounded by 

k n ~ 1+5 h{\W\) = k n - 1+s h(k\z - z'\) 

where h(t) is bounded by Ct~ n+1 for t small and is rapidly decreasing in t for t 
large. A similar argument to that for A 3 shows that the kernel of A4 is bounded by 
min(|z — z'\~ n+1 , \z — z'\~ n ~ s ), which shows that it is bounded on L p for 1 < p < 00. 

Finally, for A$, localized to k < x', we use coordinates (s — x/x',x,k' = 
k/x',y,y'). Here y — y' and s — 1, or equivalently x(z — z r ), define the diago- 
nal and a; is a boundary defining function for bfo on the support of A5. One can 
suppose that A 5 is supported in {x\z — z'\ < 1}. The kernel A 5 (k) is therefore 
bounded by (k/x)~ 1+s x n ~ 1+s h(x\z — z'\) where h is as above. Integrating in k 
from to a; (since A 5 is supported where k < x) gives a bound 



f A 5 {k)dk < x n+s h{x\z- z'\). 
Jo 



This is bounded by min(|z — z'\ n+1 , \z — z'\ ™ 5 ), which shows that it is bounded 
on L p for 1 < p < 00. This completes the proof of the proposition. □ 

5.2. Boundedness of the Riesz transform. Recall from the Introduction that 
we define the Riesz transform T by T = do (P > ) -1 / 2 . This can be written, at least 
formally, as 

(5.1) T=-d [ i?(fc)n>dfc 

Jo 
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where n> is the spectral projection onto (0, oo) for the operator P. 

Using the elementary bound ||iJ(fc)II>||i2_ > £,a < k~ 2 , we see that the integral 
converges on [ko, oo) for any ko > 0, at least as an operator L 2 (M) — > H^ 1 (M). To 
make sense of the integral near zero we use the results in this paper on the resolvent 
kernel R(k) as k — > 0. 

To compute the integral (|5 . 1 [) we divide it into pieces and also compare it to 
the classical Riesz transform. We recall that P has at most finitely many negative 
eigenvalues > — k\ > — k 2 > ■ ■ • > —k 2 M . We write Ilj for the projection onto 
the —kj eigenspace. Choose a smooth monotone function f(t) that is equal to 
for t < —kl/2 and 1 for t > -fef/4. Thus f(P) = U > . We also choose a smooth 
function ip{k) so that ip(k) = 1 for k > kl/2, ip(k) = for k < fc 2 /3. 

Consider the integral 

(5.2) -d ip(k)R(k)IL > dk = -d ip(k)R(k)f(P) dk, 
k Jo 71 Jo 

where we inserted the cutoff function i\) into integral (|5.1[) for T. Due to the 

conditions on /, the function 

9(A) = /(A) jH^fcJ—L^dfc 

is a smooth function of A which is a classical symbol of order — 1 as A — > oo. By the 
symbolic functional calculus of [5], the operator g{P) is a scattering pseudodiffer- 
ential operator of order — 1. Hence dog(P), which is equal to (|5.2p . is a scattering 
pseudodifferential operator of order 0. Using Stein [5T], chapter 6, section 5, this is 
bounded^ on D>{M) for 1 < p < oo. This treats the integral for (|5.1[) for large k. 

Now consider the operator <fiRo(k)<f>, where Ro(k) is the resolvent for the free 
Laplacian Ao on the exact cone (0, oo) x dM and is a cutoff function, equal to 
zero on {x > e} and 1 on {x < e/2} for some small e > 0. If e is small enough, this 
kernel may be regarded as living on M 2 sc and it satisfies the properties listed in 
Section ETTI (see Remark |4.2[) . We can apply the argument above to Ro(k) and we 
see that d o g(A ) is bounded on L P (M) for 1 < p < oo. 

Now consider the remaining integral, where 1 — tp(k) is inserted into the integral 
(|5.1j) for T. Choosing fco = k 2 /3 we can write this as 

2 r ko / M \ 

-d (1 - ip(k))R(k) ( Id - n i ) dk - 

It is well known via Agmon estimates pQ that the eigenfunctions corresponding to 
the Tlj, j > 1, together with all their derivatives, are exponentially decreasing at 
infinity. It follows that the Hj in the integral above contribute an operator that is 
bounded on L p for 1 < p < oo. The remainder is 

2 f k ° 

-d / (l-ip(k))R(k)dk. 
n Jo 

Rather than study this directly, we subtract off the free resolvent kernel, and con- 
sider 

(5.3) ^ J (d{R{k) - (j)R (k)(f>) + [d, (f>}R (k)^ dk. 

By Theorems ll.2l and ll.3l R(k) is in the calculus with pseudodifferential order —2 
and with index sets > at zf, > n/2 — 2 at rbo and lbo, > —2 at bfo, > at sc 
and trivial at all other boundary hypersurfaces. Similarly, our analysis shows that 



The result in Stein is stated for Euclidean space; however, by localizing in dM one reduces 
to this case. 
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[d, <j)]Ro{k)(j) and (pRo(k)<p are in the calculus with pseudo-differential order —2 and 
with similar index sets, except that [d,<f>] is compactly supported so [d, 4>]Ro(k)(p 
vanishes at all order at all boundary hypersurfaces but rb, rbo,zf. The point of 
subtracting the two resolvents is that they have the same leading term at bf a and sc, 
since these models are determined purely by the metric h(0). Hence the difference 
has index set > — 1 at bfo and > 1 at sqj- 

Now consider what happens when d is applied on the left. The operator d has 
components d Zi each of which lift from the left factor to be of the form Pbf /Oib PibPbf 
times a vector field tangent to the boundary of M| sc , and transverse to the diago- 
nal. Therefore, it increases the order of vanishing to n/2 — 1 at lbo and to at bfo. 
At rbo, by contrast, it acts globally in the z variable. So we get an improvement 
at rbo if and only if the leading coefficient is annihilated by d, or equivalently is 
constant in z. (It also increases the pseudodifferential order to —1). The lead- 
ing term of d(R(k) — 0i?o(fc)0 ~ [d, (f)\Ro(k)(f>) at rbo is at order n/2 — 2, this is 

n 2 2± 2 — 2 

c?i? r 2 bo — 4>{z)dR^ rbo where i? rb is the leading term of Ro(k) at rbo- It is impor- 

— — 2 ^—2 

tant to notice that this term vanishes if and only if R \ and i? n 2 , are constant in 

J I DO U,IL>o 

the variable z, which is exactly the case when V = and M has one end. Indeed, 
we have seen in Sections 13.71 and 1 4.51 that i? r 2 bo , as a function of z, is a multiple of 
a bounded function annihilated by P. This is uniquely determined up to constants 
when M has one end, and it is constant if and only if V = (see Remark 13. 2|) . 
When M has more than one end, on the ith component of rbo it is a multiple of 
the bounded solution of Pv = that tends to 1 at the ith end and at all other 
ends (see Remark 14. ip . In particular, in this case it is not constant. When it is 
constant, it means from Subsection 14. 51 that the next leading term at rbo is at order 
fi - 1, so we can now apply Proposition 15.11 to deduce that (|5.3|) is bounded for 
1 < p < n when V\>1 and for 1 < p < n/(n/2 — v\) := p max otherwise. Combined 
with the earlier results about the integral for large k, we have proved that T — <\>Tq§ 
is bounded on L p for 1 < p < p max , where To is the Riesz transform on the exact 
cone K + x dM. When M is asymptotically Euclidean, the cone is Euclidean R" 
and the boundedness of To for 1 < p < oo is classical. For the general conic case, 
it is shown in [lOj that To is bounded on LP for 1 < p < p max so we conclude that 
T itself is bounded on LP for the stated range. 

Finally we claim that the range of p is sharp. In this case, we can write the 
resolvent kernel near the right boundary (now as a multiple of the scattering half- 
density rather than the b-half density, which gives an extra factor of (xx') n ^ 2 ) as 

( } = g_ a «'g n /2-i(«2 ( ) + k u x -l ^xOO g/tsf-l,, {z) 

K> r(n/2 - 1)2»/ 2 - 2 oW r(z/i)2 l 'i- 1 lW 

-x r(f-i)2f- 2 ai+ r(i/ 1 )2-i- 1 vA) 

modulo O{pX 2+U2 ) + 0{p^ logp rbo ). So we have 
dR(k) = x' n ~ 2 K K~*-? dix^-'v^ (z)) 



T(f - 1)2 4 



ctually, there could be a log term at order —1 at bfo. However, this would make no difference 
to the argument below, so we ignore this detail to simplify the exposition. 
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(since d acts in the left variable z here) modulo 0(p* h ) + 0{p^ 1 logp r b )- Inte- 
grating in fc, we change variable to k! — k/x' (which gives us an extra power of x') 
and obtain 

x' n - l d{x^-\ V0 {z)) j[ _ dn' 

+ ^'*" /2 -'»*.<~'»/ h^0 dK ' 

modulo 0{x'^ +u ' 2 ) + 0(x' n logs') for the asymptotics of T at x' = 0. Since 
x n / 2 ~ x Vv (z) is constant if and only if V = and M has one end, d(x n / 2 ~ 1 ipj(z)) 
is identically zero if only and if we are in this case. In addition, the k' integrals do 
not vanish, indeed the function (for v > 0) 

F(k') := k' V K v {k') 

is smooth on (0, oo), with F(0) = 2 lJ - 1 T{v) and F(k) = F(0) + 0(n 2 ) + 0{n 2v ) as 
k — > 0, and we have the identity for v > 



/ F{n')dn' = -2 lJ - 1 ix^T{v+- 



2' 

since z u K u {z) = c u ? t ^ z ((l + t 2 )~ v -i) with c„ := -T(v + \)2 v - 1 it-^. 

Hence T = a(z,y')x' n ~ 1 + 0{x' Ul+2 ) + 0(x' n logs') at x 1 = where a does not 
vanish if and only if V ^ or M has more than one end. It follows immediately 
that the upper threshold for p is sharp in that case. A similar analysis at the left 
boundary shows that T = b(y, z')x n ~ 2+ra + 0(x n ~ 1+m ) as x — > where b does 
not vanish, showing that the lower threshold is also sharp. 

Now ifV = and M has one end, T = c(z, y')x' Ux+ % +0(x' U2+ %)+0(x' n ~ l logs') 
with c{z,y') = d z (x^~ 1 v^ 2 (z,y')) and v V2 solving P(x^~ 1 v U2 (z,y')) = and 
v V2 (x,y,y') ~ 7(2/', y)x~ 1 ' 2 as s — > for some 7(2/,?/') 7^ 0. Thus c ^ and 
this proves the upper range is sharp in that case, at least when v\ < i/q + 1 = ^. 

This completes the proof of the theorem. 



6. Appendix: proof of Proposition 12. 1(J1 

This result will follow from the existence of a 'triple space' M| sc having the 
property that there are b-fibrations down to M 2 sc which are lifts of the three 
projections M 3 x [0,k ] -> M 2 x [0,fc ]. We now proceed to define such a space. 

We first define M| b to be the 'total blowup' of the space M 3 x [0, fco] in the 
sense of [7], that is, with the codimension 4 corner blown up, followed by the 4 
codimension 3 corners (which may be done in any order as they are separated after 
the first blowup), followed by the 6 codimension 2 corners (which likewise may be 
done in any order). Then the three projections from M 3 x [0, fco] M 2 x [0, fco] 
extend from the interior to b-fibrations from b — » M 2 b (this can be deduced 
from the lemmas in section 2.4 of [7], for example). Let us denote these b-fibrations 
7r^» for * = L,C,R (left, centre and right) according as it eliminates the right, 
centre or left variable respectively. 

The space M 2 sc is the space M| b with the boundary dAk,b at bf blown up. 
The inverse image of this submanifold in Mj 3 b under any one b-fibration 717,,* is the 
union of 2 p-submanifolds, J* and G* . Their intersection properties are such that 
once the intersection K of the G* is blown up, the lifts of G* are all separated, 
while the J* only intersect the corresponding G*. Therefore, we define 

(6.1) M ksc = [M k y, K-. Gl,Gc, Gr; Jl, Jc, Jr] ■ 
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Lemma 6.1. The three projections from the interior of M 3 x [0, ko] to the interior 
of M 2 x [0, ko] extend by continuity to b-fibrations from M 3 sc to M% sc . 

The proof uses the following properties of blowups and b-fibrations (see [H] or 
®)- 

Lemma 6.2. Suppose that X and Y are manifolds with corners and f : X —>Y is 
a b-fibration. Let S andT be p-submanifolds ofX. If either 

(i) ScT, or 

(ii) S is transverse to T, 
then [X;S;T] = [X;T;S]. 

(Hi) Suppose that S is b-transversal to f and such that f(S) is not contained in 
any codimension 2 face of Y . Then f restricted to X\S extends to a b-fibration 
from[X;S]^Y. 

Proof of Lemma \6.1\ The proof runs along the same lines as the proof of the ex- 
istence of b-fibrations from the triple scattering space M 3 C to M 2 C in Section 23 of 
[IB] , and we follow this proof to the extent possible. We refer to [17 for definitions 
of terms such as b-fibration and b-transversal. 

Since the spaces M| b and Mf sc are symmetric under permutation of the M 
factors, it is enough to show that one of the projections from the interior of M 3 x 
[0, ko] to the interior of M 2 x [0, ko] extends to be a b-fibration. We shall take the 
projection ttr that projects off the right factor of M. 

We first show consider the space M| b . We temporarily use 4-digit binary codes 
for the faces of M 3 x [0, ko] where 0000 represents (dM) 3 x {0}, 1010 stands for 
M x dM x M x {0}, 0011 stands for (dM) 2 x M x [0,fe ], etc. Then 

M| b = [M 3 x [0, ko]; 0000; 1000, 0100, 0010, 0001; 1100, 1010, 1001, 0110, 0101, 0011]. 

Here the four codimension 3 faces may be permuted among themselves since they 
become disjoint after the 0000 blowup, and similarly the six codimension 2 faces 
may be permuted among themselves. So we can write 

M| b = [M 3 x [0, fc ]; 0000; 0010, 0001, 0100, 1000; 1010, 0110, 0011, 1100, 1001, 0101]. 

Now the 0000 and 0010 blowups may be commuted using (i) of Lemma \6. 21 Then 
the 1010 blowup may be moved to the left past 0001 (disjoint, once 0000 has been 
blown up), 0100 (disjoint), 1000 (contains) and 0000 (contains). Similarly 0110 and 
0011 may be moved to the left, obtaining 

M| )6 = [M 3 x [0, fc ]; 0010; 1010, 0110, 0011; 0000; 0001, 0100, 1000; 1100, 1001, 0101]. 
The first four blowups here give us M| b x M. So we have shown that 
Ml b = [Af£ 6 x M\ 0000; 0001, 0100, 1000; 1100, 1001, 010l]. 

= [Af£ 6 x M; bfo x dM- bf x dM, lb x dM, rb x dM; zf x dM, rb x dM, lb x dM] . 

Now we include the blowups that turn M 3 b into Af s 3 c b . Note that these may be 
reordered Gr, Jr, K, Gc, Jg, Gl, Jl since K C Gr and K and Jr are disjoint. So 
we have 

M Lc = [ M l.b x M; bf x dM; bf x dM, lb x dM, rb x dM; 

zf x dM, rb x dM, lb x dM; G R ; J R ; K; G c , J C \ G L , J L ] ■ 

Recall that Gr is the submanifold dhi^k,b X dM C M 2 b x M, and Jr is the 
submanifold dhi&k,b x M C M 2 b x M, lifted via the blowups listed to their left. 
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Since dbfAk,b C M| sc does not intersect zf, lbo, rbo, lb, rb, they can be commuted 
to the left past five blowups. So we get 

Ml sc = [Ml b x M ; bf x dM; bf x dM; G R ; J R ; lb x DM, rb x 0M; 

zf x 0M, rb x dM, lb x 0M; if; G c , J C \G L , J L ] . 

We can commute Gr to the left of bf x dM since it is contained in bf x dM; then 
we can commute Jr to the left of bf x dM since these two submanifolds become 
disjoint after Gr is blown up; and finally we can commute Jr past Gr since it 
contains Gr. Finally, Jr and bfo x dM are transverse so these two blowups can be 
commuted, using (ii) of Lemma lG.fl We end up with 

M| iSC = [Ml b x M ; bf x dM; J R ; G R ; bf x dM; lb x dM, rb x dM; 
zf x 0M, rb x dM, lb x 0M; K;G Cl Jc]G L , J L ] 
= [Ml iBC x M; bf x 0M; G fl ; bf x 0M; lb x 0M, rb x 0M; 

zf x dM, rb x dM, lb x 0M; X; G c , J C \ G L , J L ] ■ 
Now using (iii) of Lemma 16.21 repeatedly we see that there is a b-fibration from 
M fe, sc t0 M l,sc- This completes the proof. □ 

We resume the proof of Proposition 12.101 We only prove the result in the case 
that (a sc , <2bf j a zi) — (0,0,0) since the general case follows easily from this. Com- 
position may be defined in terms of pullbacks and pushforwards. We first need to 
take care of the density bundles. We have an isomorphism 

(6.2) n* L nl /2 (Ml sc ) ® n* c & b /2 (Ml sc ) ® 7r^ /2 (Aff !SC ) = fi b (M fc 3 sc ) ® \f | 1/2 

where &b(Ml ) is the lift of f2;,(M| 6 to M? sc . This allows us to express the 
composition A o B as 

(A o = (7r c )„ (tt^A (g) n* R B <g 77^.1/ ® I — p 1/2 ^j 

where ^ is a b-half-density on M^ sc . 

Given A and B as in the proposition, decompose them into A = A& + A p h g , 
B = £>a + ^phg where Aa and B& are supported close to A/j jSC and classical 
conormal to it, uniformly to the boundary, while j4 p h g and -B p h g are smooth in 
the interior and polyhomogeneous at the boundary with the given index sets. The 
composition result for A p h g and -B p h g follows directly from the b-fibration property 
of Lemma E2] and the Pushforward theorem (Theorem 5 of [I?]). For example, the 
second line of (|2.14p follows since there are two boundary hypersurfaces of M| sc 
that project to zf under ttq, one of which projects to zf under both 7Tl and 7TR, and 
the other of which projects to rbo under tt^ and lbo under ttr. 

Next, the result for A& and Ba follows from the geometry of the lifts of Ak, S c C 
M| sc to Ml sc by each of the 3 b-fibrations. Namely, each of the lifts is an interior 
p-submanifold, each pair of which intersect transversally at the lifted triple diagonal 
which projects diffeomorphically to Afc iSC under ire- 

Finally, consider the composition of Aa with S p h g (the case j4 p h g with Ba is 
equivalent). We observe that the composition property holds if A is a differential 
operator of any order. We can write A a = A' o D + A" , where D is an differential 
operator of order 27V + m which is elliptic as a conormal distribution to A^, A' 
is supported close to Afe iSC and is conormal of very negative order —N, and A" is 
smooth and supported close to Afc iSC . The composition of A" with -B p h g can be 
treated using the Pushforward theorem. The differential operator D maps S p h g to 



^This step is the only essential addition to the arguments of 1161 
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another polyhomogeneous kernel with the same index family. On the other hand, 
A' has a kernel which is C N ~ n ~ 1 . Therefore A' o D o B p h g is polyhomogeneous to 
some finite but large order ~ N. Since N is arbitrary, this shows poly homogeneity 
of A A o B phg . 
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